A measurement scheme for purity based on two two-body gates 
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A scheme for measuring the purity of a quantum system with a finite number of levels is presented. 
The method makes use of two VSWAP gates and only hinges on measurements performed on a 
reference system, prepared in a certain pure state and coupled with the target system. Neither 
tomographic methods, with the complete reconstruction of the state, nor interferometric setups is 
needed. 
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I. INTRODUCTION 

Purity is a good measure of the coherence of a quan- 
tum system. It is unity for a system in a pure state 
(complete quantum coherence), while it reduces to 1/N 
for an A^-level system in a completely mixed state (no 
quantum coherence). Purity has played a central role in 
the discussion of the quantum measurement problem [l[ , 
in the sense that the occurrence of decoherence, consid- 
ered to be responsible for the transition from a quantum 
to a "classical" (mixed) state by the action of measure- 
ment, is one of the main issues to be resolved within the 
framework of quantum theory 0. Furthermore, its im- 
portance is well recognized in the field of quantum infor- 
mation, communication and computation, for the purity 
of a subsystem is an important measure of the entan- 
glement of the total bipartite system 0] : the higher the 
entanglement between the two subsystems, the lower the 
purity of the reduced density matrix of the subsystems, 
and vice versa. 

Although purity is an important physical notion, it 
is not so simple to devise an experimental procedure to 
directly measure it. Remember that the quantum me- 
chanical expectation value of an operator O in the state 
p is expressed as Tr{Op}. This is a linear functional of 
p, while purity II = Ti{p 2 } is a quadratic functional. Pu- 
rity is usually calculated only after the state p has been 
reconstructed by some tomographic methods There 
are, however, interesting proposals @-0| to measure lin- 
ear and nonlinear functionals of the density matrix, in- 
cluding purity Tr{p 2 }, without resorting to any tomo- 
graphic methods. One assumes that copies of the quan- 
tum system under scrutiny (target system) are available, 
so that the state p <£> p <E> ■ ■ ■ can be fed to an interfer- 
ometric setup in a quantum network. In particular, to 
determine the purity Tr{p 2 } of the target system, one 
requires a controlled-SWAP operation and an ancillary 
two-level quantum system (qubit). 

In this paper we present an alternative scheme for mea- 
suring the purity of a quantum system with a finite num- 
ber of levels [typically a two-level system (qubit) or a 
three- level system (qutrit)], without resorting either to 



the state tomography of the density operator or to inter- 
ferometry in quantum networks. In this way, a relatively 
small number of (different) measurements arc required 
Q. In particular, our scheme consists of two simple two- 
input gates, such as VSWAP, and does not rely on the 
three-input controlled-SWAP gate [HIH, whose explicit 
construction is not trivial. Moreover, the necessary ancil- 
las, that interact with the target systems [H, Q, could be 
automatically supplied by projective measurements per- 
formed on a fraction of the ensemble itself of target sys- 
tems. As we will show in explicit examples, our strategy 
does not rely on the independent preparation of ancil- 
lary systems, provided the ancillas are still available for 
further manipulations after the projective measurements. 

This paper is organized as follows. In Sec. [TTl we 
present the general ideas. We show here how to extract 
information about the purity of a quantum system. As 
stated above, since purity is a quadratic functional of the 
density matrix, we prepare the initial state in a tensor- 
product form which includes two system density matri- 
ces. An ancilla system is introduced and coupled to the 
system under consideration. A unitary operation (gate), 
coupling the system and the ancilla, is applied twice in 
succession. Finally, a measurement on the ancilla yields 
information on the purity of the system. Neither tomo- 
graphic nor interferometric setups is introduced. These 
general ideas are then made concrete by explicitly con- 
structing the unitary operator and by identifying the nec- 
essary measurements in the case of a two-level system 
(qubit) in Sec. a three-level system (qutrit) in Sec. 
ITVl and a general AT-level system in Sec. |V] Section [Vl] 
is devoted to the summary and discussion. An Appendix 
is added to describe, in terms of a "generalized" Bloch 
vector, how the scheme works in general. 



II. GENERAL IDEA AND FRAMEWORK 

Let a quantum-mechanical system (target system) 
with N discrete levels be described by a density matrix 
p, which can be characterized by a generalized (TV 2 — 1)- 
dimcnsional Bloch vector a. Our goal is to measure its 
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purity II = Tr{p 2 }. Assume that we can prepare the tar- 
get state p in duplication, i.e., p® p. We also prepare 
another quantum system (ancilla) with N discrete levels 
in a particular pure state w. The initial state reads 



• p <x> oj. 



(1) 



Observe that this ancilla can be supplied by a fraction of 
the ensemble of the target state p by performing an ap- 
propriate projective measurement, which is required in 
our scheme for the estimation of purity, as we will sec 
below, and one could thus reduce the number of different 
resources required. We hrst make one of the duplicated 
target systems interact with the ancilla via a unitary evo- 
lution operator (gate) U. After the interaction, the infor- 
mation about the target state p has been transferred to 
the ancilla and they are in an entangled state. The total 
system is then exposed to another interaction, this time 
only between the other target system and the ancilla, 
governed by the same type of unitary evolution operator 
U. Finally, a suitable observable of the ancilla is mea- 
sured, from which information about the purity of the 
target system can be extracted. 

The state of the ancilla just after each step, and 
w( 2 >, is expressed in the following way. Define the map 
acting on two quNits 



A(p 1 ®p 2 )=Tr 1 {U{p 1 ®p 2 )U^}, 



(2) 



where U is a unitary operator on C N ® C N and the trace 
Tri is taken over the first system (left-hand side of the 
tensor product). The states of the ancilla after the first 
and second step read 



and 



= A(p® cj) 



oj {2) = A(p® A(p®u)), 



(3) 



(4) 



respectively. Note that we do not measure the target 
systems (and therefore take the partial traces to get the 
ancilla state). 

Since the final state of the ancilla, represented by the 
reduced density matrix oj^ 2 \ depends quadratically on 
p and therefore quadratically on its Bloch vector a, the 
expectation value of an ancilla's observable is expected 
to carry part of the information about a 2 (because in 
general there remains no larger symmetry that keeps a 2 
invariant). This is why we have presented the target 
state p in duplication. The purity of the target state p 
can now be evaluated, eventually by supplementing addi- 
tional information about vector a, that can be obtained 
by ordinary projective measurements directly on the tar- 
get state p. A schematic diagram is shown in Fig. [1] 

To be more definite, assume that the ancilla is prepared 
in a pure state 
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FIG. 1: Schematic diagram of the procedure. Evolution due 
to Eqs. Q and @ : the output of the first gate, whose inputs 
are p and oj, is fed to the second gate, together with the other 
p. A projective measurement O is performed on the ancilla, 
but not on the target states (that are traced away). 



and we measure the probability of finding the ancilla in 
the same state 



O = 1 ® 1 ® \n)(n\ 



(6) 



Clearly, this observation is performed on the final state 
a/ 2 - 1 . Furthermore, suppose that the interaction Hamil- 
tonian H between each target system and the ancilla is 
a SWAP operator 



H(\i) ® |n)) oc \ri) ® \i), Mi,n. 



(7) 



Then it is easy to sec that, since the successive gate op- 
erations induce on the basis vector the transformation 
n) = \i) ® \ j) ® \n) — s- \n,j,i) — > \n,i,j), the expec- 
tation value of O contains terms like 

(O) ~ ^{n,i,j\(p® p®\n){n\)\i,j,n) 

= ^( n W)( i \p\ n ) 



(n|p 2 |n), 



(8) 



which are nothing but the diagonal components of p 2 . 
By collecting these results for a complete orthonormal 
set we obtain the purity 



n 



n\p 2 \n) 



(9) 



OJ 



\n)(n\ 



(5) 



This clarifies which elements are important for the above 
general framework to work properly: i) a SWAP Hamil- 
tonian, which brings about a two-body gate, and ii) a set 
of projective measurements on the ancilla. As we will sec 
below, we shall also need iii) a direct measurement on p 
to single out (n\p 2 \n) from (O). We emphasize that the 
scheme requires only two simple two-body interactions 
(SWAPs) and does not rely on more complicated gates 
like the controlled-SWAP. 

The general scheme described above will now be tested 
on two simple examples: a two-level system (qubit) and 
a three- level system (qutrit). The extension to general 
iV-level (./V > 3) systems is straightforward and will be 
presented in Sec. |Vj Incidentally, as is clear from the 
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above discussion, in practice one needs to calculate the 
reduced density matrices at each step as functions of the 
Bloch vector a, or alternatively, one can directly compute 
the expectation value (O) in the final state, as is shown 
in Sec. [VJ 

III. QUBIT (N = 2) CASE 

Let a quantum-mechanical two-level system (target 
qubit) be described by a density matrix p, characterized 
by an unknown Bloch vector a, 

p =±(1 + a- a), (10) 

where <x is the vector of Pauli matrices. The square mod- 
ulus of a is related to the purity II of p by 

n = TiV} = ^(i + « 2 ), ^< n <!- (ii) 

If the system is in a pure state, II = 1 and a 2 = 1, 
while for the completely mixed state p = 1/2 purity is 
minimal II = 1/2 and a = 0. We introduce an ancilla 
qubit, prepared in a pure state 

u = \n){n\ = ^(l + n-a), (12) 
(n-a)\n) = \n), \n\ = 1 (13) 

and make it interact with the target qubit. Assume that 
the unitary operator (gate) U representing the interac- 
tion between each of the target qubits and the ancilla is 
given by 

U= 4=(1L® 1-iS), (14) 
v2 

which can be generated by the SWAP Hamiltonian 

1 1 3 

S = -1(8) 1 + -£V<8><t* (15) 

with the Heisenberg interaction as U = e — 'ff* 5 | fff=7r / 4 . 
Indeed, by noting the property of the SWAP operator 
S 2 = 1 <g) 1, one gets U(t) = e~ lgts = cos(gt)l <g) 1 - 
iSsin(gt), which reduces to (fT4"]) at t = ir/Ag. Observe 
incidentally that this unitary operator is just a realization 
of vSWAP, since its square yields U 2 = —iS, which is 
nothing but the SWAP 0. 

It is easy to see that if pi = i(l + a ■ cr) and p 2 = 

i(l + 6-<r), the map p 2 1] = A(p x ® p 2 ) = ±(l + c-<r) is 
mirrored into the map of the corresponding Bloch vectors 

c = A(a,b) = i(a + b + a x b). (16) 

In the scheme outlined in Fig. [TJ pi and p 2 are given by 
(JTOJ), while LiJ is given by (fT2")) . so that the Bloch vector of 



the state of the ancilla after the first step a/ 1 ' = ^(1 + 
bi ■ cr) reads 

b\ = A(a, n) = i(a + n + a x n). (17) 

From @, PB]) and (JT7J) , the Bloch vector b 2 character- 
izing the reduced density matrix of the ancilla after the 
second interaction a/ 2 ) = i(l + 62 • cr) is given by 

6 2 = A(a, 61) = A(a, A(a, n)) 

= ^[3a + n + 2axn + ax (ax n)]. (18) 

Observe that this expression is quadratic in the Bloch 
vector a. If we measure the spin of the ancilla along 
direction n in the final state ui^ 2 \ its expectation value 
N is given by 

Af = Tr{(n-o-)uj i2) } = n-b 2 

= j[l + 3n • a+ (n ■ a) 2 -a 2 }. (19) 

The purity of the target state p can thus be written as 
3 1 

n = 1 + -n • a + -(n • a) 2 - 27V . (20) 

Evaluation of this expression requires knowledge of n ■ a, 
which can be easily obtained by directly measuring the 
spin of the target system p along n, 

Tr{(n • cr)p} = n ■ a, (21) 

or by measuring the spin of the ancilla a/ 1 ) along n, 

Tr{(n • er)w (1) } = n ■ b x = -(1 + n ■ a). (22) 

Both these measurements can be performed indepen- 
dently, e.g. by using a portion of the targets or ancillas. 

Notice that no tomographic method involving the com- 
plete reconstruction of the state p has been invoked; only 
a single type of measurement on the ancilla, that is, the 
measurement of the spin along direction n, is required to 
obtain the purity of the target state p. In other words, 
there is no need to collect data with different n's. Notice 
also that n can be chosen arbitrarily. 

Furthermore no intcrferometric setup in a quantum 
network d, @ is required, significantly reducing in this 
way the necessary experimental steps. Indeed, only two 
simple VSWAP gates are employed. 

Finally, the ancilla could even be supplied from the 
ensemble of the target state p, when one performs the 
measurement on a part of it and finds it in the same 
direction n, since the state after the measurement with 
an affirmative result is reduced to the state w = |n)(n|. 
This also yields the necessary information, i.e. the quan- 
tity n ■ a, related to the probability of finding the target 
system in state u>. 
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IV. QUTRIT (N = 3) CASE 



and can be shown to be characterized by the Bloch vector 



The extension of the above scheme to multi-level sys- 
tems is straightforward (although not trivial). In this 
section, we consider N — 3 case (qutrit). Generaliza- 
tion to still higher N is given in the Appendix, while a 
different derivation is given in the following section. 

The density matrix of a three-level system can be writ- 
ten as 



1 



-1 



1 

71 



a c X c 



1 



■1 



1 



v/3 



a ■ A, 



(23) 



where, here and in the following, the summation over 
repeated indices (in this case from 1 to 8) is implicit, 
a is a generalized eight-dimensional Bloch vector and 
the generators of the SU(3) group X a /2 (a = 1, ... ,8), 
that are all traceless Tr{A a } = 0, satisfy the commu- 
tation and anti-commutation relations with the totally 
anti-symmetric and symmetric structure constants 

[A a ,A b ] = 2if abc X c 1 {A a ,A b } = -6 ab + 2d abc X c , (24) 

3 

with the normalization condition 

Tr(A a A b ) = 2S ab . 
The purity of the system is 

H = Tr{p 2 } = \ + \a 2 

and a pure state with purity II = 1 is again characterized 
by a Bloch vector of unit length \a\ = 1. We have to pay 
due attention to the fact that not all the states expressed 
as in (|23jl arc physically acceptable. In fact, in order to 
insure the positivity of the density matrix, the physical 
domain of the Bloch vector is much more restrictive than 
\a\ < 1 Q. For the moment, we assume that only physi- 
cally acceptable Bloch vectors have been chosen when we 
write the density matrix in terms of the generators A a /2. 
We shall return to this point later. 

We now let a target qutrit, whose density matrix p is 
characterized by an unknown Bloch vector a as in (|23|) . 
interact with an ancilla qutrit, prepared in a prescribed 
pure state 



(25) 



(26) 



f 



-1 



f 



vr 



i. 



(27) 



Assume that the interaction is again vSWAP as in (|T4| , 
with the SWAP Hamiltonian S for a couple of qutrits 

m 



s 



1 



-1 



i + l -x< 



A< 



(28) 



(Additional details on this general structure can be found 
in the Appendix.) After the unitary evolution engen- 
dered by this operator, the reduced density matrix of the 
ancilla reads 



1 



1 



(29) 



6i = A(a, n) = i (a + n + x n ) 



where a simplified notation has been introduced for the 
anti-symmetric product between eight-dimensional vec- 
tors 



(a x n) a = f abc a b n c , ax n = -nx a. (31) 

Finally, the reduced density matrix oj^ of the ancilla 
after the second interaction is characterized by the Bloch 
vector 62 

b 2 = A(o, bi) = A(a, A(o,n)) 

13 1 1 

= —n H — a H —(a x n) H — a x (a x n). (32) 

4 4 ^ y 1 3 

By measuring the operator n ■ X on the ancilla, we get 
the expectation value 

N = Tr{(n • A)w (2) } = ^=n • b 2 
V3 

= ^ + \{n.a)-l { axnf). (33) 

Notice that the last term contains a bilinear factor in 
a, which can be rewritten as 

(a x n) 2 

= f abc f cde a a n b a d n e 
2 

= - [a 2 - (n ■ a) 2 ] - (n * a) 2 + (a* a) ■ (n* n), (34) 
o 

where a symmetric product between two eight- 
dimensional vectors has been introduced, yielding an- 
other vector 

(a*n) a = d abc a b n c , 

a*n = n*a, n ■ (a * n) = a ■ (n * n), etc. (35) 

Recall that in dimension D (= N 2 — 1) higher than three, 
the magnitude of a D-dimcnsional vector a cannot be 
simply determined by its component a ■ n along a partic- 
ular unit vector n, and by one of its normal components 
I a x n|, since there are more than one directions normal 
to a plane spanned by two vectors in D dimensions, which 
can alternatively be seen as the breakdown of the vector 
relation n x (a x n) = a — (a ■ n)n, valid in three dimen- 
sions. This means that the above expectation value (|33[) . 
supplemented with the value of n ■ a, is not enough to 
extract the Bloch vector squared a 2 . One may overcome 
this problem by measuring not a single particular opera- 
tor n- A fixed by a given n, but several (in practice three) 
operators specified by properly chosen unit vectors. At 
this point, we have to be careful in choosing n so that 
the resulting density matrix really represents a physical 
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state [9(. In order to trivially satisfy this physical-state 
condition, we choose the following unit Bloch vectors 

m = (0,0, V3/2, 0,0, 0,0, 1/2), 

n 2 = (0,0,-V3/2,0,0,0,0,l/2), 

n 3 = (0,0,0,0,0,0,0,-1), (36) 

which correspond to the states of the ancilla represented 
by the diagonal density matrices 

1 0\ /0 0\ /0 0\ 
000 , 010 , 000 , (37) 

oooy \oooy \o o iy 

respectively Choosing one of these unit vectors simply 
means that we prepare the ancilla qutrit in one of the 
three levels and measure its population. We easily evalu- 
ate the quadratic term in a by inserting the actual values 
of the structure constants 

(axn 1 ) 2 = |((« 1 ) 2 + (« 2 ) 2 + (« 4 ) 2 + (« 5 ) 2 ), 

(a x n 3 ) 2 = \ ((a 4 ) 2 + (a 5 ) 2 + (a 6 ) 2 + (a 7 ) 2 ). (38) 

It is evident that these values, together with the values 
of a 3 and a 8 , that are also available just by measuring 
the populations of the three levels in p (or in w^), are 
enough to determine the Bloch vector squared a 2 . In 
other words, the measurements of the populations of the 
three levels (in states u/ 2 ) and p or a/ 1 )) are enough to 
determine the purity of the qutrit system p characterized 
by the (eight-dimensional) Bloch vector a. 

V. GENERAL N CASE 

The qutrit case of the preceding section can be 
straightforwardly generalized to higher N . however, the 
results are a bit involved and are collected in the Ap- 
pendix. In this section we follow a simpler route and 
obtain the purity of a general iV-level system, according 
to the general idea presented in Sec. [TTJ 

Recall that the density matrix p can be expressed as 

P = ^2p a \tp a )(ip a \, ^2p<* = TrM = 1) (39) 

a a 

in terms of its eigenvectors \tp a ) belonging to the eigen- 
values p a (0 < p a < 1), which are orthonormal to each 
other (ipalipp) = <W- We prepare the ancilla state oj 
in one of the iV-levels, say \n) (n = 1, . . . ,N), so that 
lu = \n)(n\. Since our initial state is of the form (fTJ), we 
first calculate the action of two successive vSWAP op- 
erations (fT4|) on the state \i/j a ) <8> \i>p) ® \n) = \a,P,n), 
obtaining 

\a,P,n) -> i[|a,/3, n) - i\a,n,0) - i\n,P,a) - \n,a,P)\, 

(40) 



where the VSWAP operation on the first and last (an- 
cilla) entries is followed by that between the second and 
the last (ancilla) entries. If the projective measurement 
on the state \n) is performed on the ancilla (i.e., the last 
entry), this state vector reduces to 

h\a, (i) - i(n\P) \a, n) - i(n\a)\n, (3) ~ (n\/3)\n, a)}. (41) 

Thus we know that the probability of finding the ancilla 
in state \n) in the final state after the two successive 
VSWAP operations on the initial state p®p®to is simply 
given by 

\^VaPfj^{[\a, (i) -i(n\0)\a,n) 

- i(n\a)\n, (3) - (n\P)\n,a)] 
x [{a,P\ +i(n\P)*(a,n\ 

+ i(n\a)*(n,p\-(n\p)*(n,a\]} 

a 

+ £$>«K»hMl a + 5(£p«IW«>l 3 ) a 

a a 

= \-\^{p 2 \n){n\} 

+ l^{p\ n ){n\} + \(Tr{p\n){n\}) 2 . (42) 

Therefore, if one sums up the result for each measure- 
ment, from n = 1 to JV, one ends up with 

JV 

^-^TrVj + i^^pIn)^}) 2 . (43) 

" n=l 

[In the case N = 2 (qubit case), one easily finds that 
Tr{p 2 |n)(n|} contains the Bloch vector squared a 2 for an 
arbitrary |n) and no further measurements are necessary 
to obtain the purity] Purity, represented by the second 
term in the above expression, is thus extracted from the 
result of the measurements of the ancilla state after the 
two VSWAP operations, if it is supplemented with the 
information on the population of each level Tr{p|n)(n}, 
which is obtained by the usual projective measurement 
on p. 

VI. SUMMARY AND DISCUSSION 

The above examples explicitly show how the general 
idea presented in Sec. [II] can actually be implemented in 
practice. The scheme presented enables one to measure 
the purity of a quantum mechanical system by resort- 
ing neither to tomography nor to interferometry. This 
implies that one would need a relatively small number 
of experimental manipulations. Indeed, in this frame- 
work, one only needs to measure the population of each 
discrete level: the number of different types of measure- 
ment scales linearly with N, a situation to be contrasted 
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with that encountered in the strategy that makes use of 
state tomography, where N 2 — 1 independent elements of 
the density matrix have to be determined. 

As explained in Sec. [ill one needs to prepare the sys- 
tem in duplication in order to estimate its purity, since 
purity is a quadratic functional of the density matrix. 
It is known and is easily confirmed that the expectation 
value of the SWAP operator S, which is a two-body oper- 
ator, on such a duplicated state directly yields the purity 
of the state, Tr{S(p ® p)} = Tr{p 2 }. It is, however, 
not trivial to realize such multi-body measurements [l2| . 
If one is allowed to measure only local observables, not 
multi-body operators, one would need one more system, 
such as an ancilla, to get information on the purity. The 
situation is exactly what happens in our case and also in 
those cases where an interferometric setup is used [H, @ . 

Since we need N types of projective measurements for 
the target system and for the ancilla after the successive 
gate operations, and since the system has to be prepared 
in duplication, the number of necessary resources is 3iV, 
which is to be contrasted with the three resources nec- 
essary in the interferometric setup 0, Q . (The three re- 
sources are the target system in duplication + one ancilla: 
the "path" degree of freedom in the schemes presented 
in [1, p is actually made up of two levels of the ancilla 
qubit.) In spite of this apparently less appealing feature, 
the present scheme is free from potential difficulties that 
are inherent in the schemes based on the controlled gates 

ra. 

The ideas put forward i n M , Q were implemented in 
two experiments. In Ref. [12[, a SWAP (flip) operator 
was directly measured to obtain the overlap of photonic 
polarization states, without introducing an interferomet- 
ric setup, while in [13j the fidelity of two NMR qubits was 
measured. The scheme discussed in this article, on the 
other hand, makes use of two-body unitary gates, that 
can be easily realized by means of a familiar interaction 
Hamiltonian. Moreover, no controlled gate is necessary. 
The projective measurement at each level, required in 
the present scheme, is nothing but a measurement of the 
population of that level. 

Finally, the examples presented in this paper depend 
on the practical realization of SWAP Hamiltonians and 
this may not be a trivial problem for general A^-level 
systems, though the interaction reduces to the familiar 
spin-exchange one in the qubit (N = 2) case. In this re- 
spect, it is an open problem, yet to be clarified, whether 
the present scheme could be generalized to other Hamil- 
tonians. On the other hand, the general ideas presented 
in Sec. |TT] could be extended to more general nonlinear 
situations, e.g. in order to estimate Tr{p n } with n > 2. 
It would be worth exploring such a possibility along the 
same line of thought, since the scheme only requires sim- 
ple two-body gate operations and projective measure- 
ments and is thus relatively easy to implement in ex- 
periments. 
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Appendix: General N case: Another exposition 

We rederive here the result obtained in Scc.[V]by using 
an alternative method, with the hope that the different 
mathematical techniques employed will aid in future ex- 
ploration of the field. The derivation is more lenghty, but 
can lead to generalization for higher-order functionals of 
the density operator. 

The density matrix of an N- level system can be written 

as 



1 N-l „ „ 1 N-l rK . 

p=—l+\ a c T c =— 1 + \ a-T, (A.l) 

H N V 2N N V 2JV v ' 

where the index c runs from 1 to N 2 — 1, a is an (A^ 2 — 
l)-dimcnsional Bloch vector, and the generators of the 
SU(A) group T a /2, that are traceless Tr{T a } = 0, satisfy 

[T a , T b ] = 2if abc T c , {T a , T b } = —5 ab + 2d abc T c 

(A.2) 

with the normalization condition 

Tr(T a T b ) = 25 ab . (A.3) 

Here f abc and d abc are the SU(AT) totally anti-symmetric 
and symmetric structure constants. One can also derive 
the relations 

T a T b = -^S ab + (d abc + if abc )T c , (A.4) 



Tt j m = --5 ij 6 ke + 25 u S jk . 



(A.5) 



The structure constants are written in terms of the gen- 
erators as 

d abc = 1 Tl (r a {T b ,T c }), f abc = — Tr(T Q [T b ,T c ]), 
4 Ai 



from which the following relations follow 



(A.6) 



jabc jabd pj^cd ^abc^abd 



N 2 -4 



N 



S cd , d a 



0, 



jabc jcde ^aa^be ^bd^ae^ _j_ ^adc^ceb ^bde^ 



dc jceb jbdc jcea 

(A.7) 
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The purity of the system in terms of the Bloch vector a 
reads 



II = Tr{p 2 } 



1 

N 



N - 1 
1 

N 



(A.. 



and thus a pure state with purity II = 1 is characterized 
by a Bloch vector of unit length a 2 = 1. 

A quantum iV-level system, described by the density 
matrix (jA.ll) . with a unknown, is coupled with another 
quantum N- level system (ancilla), prepared in a pure 
state 



1 

N~ 



N-l 
2N 



-n ■ T, 



(A.9) 



by the VSWAP gate U as in dEJ) with the SWAP Hamil- 
tonian S for a couple of AT- level systems 



(A.10) 



After this unitary gate, the Bloch vector bi characteriz- 
ing the reduced density matrix ui^ of the ancilla is 



bi = -(a + n) 



N - 1 
2N 



(A.11) 



where the same notation as in the text has been 
used for the anti-symmetric product between (N 2 — 1)- 
dimcnsional vectors (a x n) a = f abc a b n c . After the sec- 
ond gate, again given by the unitary operator (fT4| . act- 
ing on the state p of another target system and the state 
u/ 1 ) of the ancilla characterized by bi (see Fig. [1]), the 
Bloch vector characterizing the reduced density matrix 
a/ 2 -* reads 



These unit vectors can certainly be written as linear com- 
binations of n3, rig, . . . , rinv-x)2_i, and h N 2_i, where 
n a (a = k 2 — 1; k = 2, . . . , N) is the unit vector point- 
ing in the direction of the (k — l)th Casimir operator, 
i.e., T a = fi a ■ T. Choosing one of such unit vec- 
tors, corresponding to one of the diagonal density ma- 
trices in (|A.14[) . means to prepare the ancilla in a pure 
state in which only one of the N levels is populated 
and to measure its population after the successive gate 
operations. In this case, the term quadratic in a in 
(IA.13|) contains terms like (a x n a ) ■ (a x hp), where 
a, p e {3, 8, . . . , TV 2 - i}. Observe that 

(axh a )iaxh p ) = i(o.T)y(a.T) ii (T5-r?.)(^-2^), 

(A.15) 

where the summations over i and j arc understood and 
the fact that both Casimir operators T a and T@ are di- 
agonal matrices has been used. 

We now show that the summation of (|A.15[) yields a 
desired result, eliminating unwanted terms like the last 
two terms on the right hand side of (f34]> , and leaving only 

relevant terms. Let n, (i = 1 JV) be the unit vector 

corresponding to the ith-level pure state whose density 
matrix is of the above diagonal form (|A.14[) with the only 
matrix clement 1 at the component. It is expressed 
as a linear combination of h a 



la '"C 



(A.16) 



and the density matrix reads 



b 2 = A(o,6i) 
1 3 



A(a, A(a,n)) 



-n + -a 

4 4 



N-l_ 
2N ' 



N 



2N 



-a x (n x a). 

(A.12) 



We measure the operator n ■ T of the ancilla in its final 
state a/ 2 ) and get its expectation value 



, w , 1 N-l 1 

1 A 1 N V 2N N 



N-l 
2N 



(i)rpC 



(A.17) 

Here the repeated Greek indices mean the summation 
over the diagonal (Casimir) operators a € {3, . . . ,A 2 — 1}. 

Since the coefficient 7^ is given by 



# = Tr{(n-Ty 2) } 



(2)1 _ 



2(N- 1) 



N 



n ■ bi 



2(N - 1) 



N 



N - 1 
2N 



(a x n) 



(A.13) 



In order to extract a 2 = ^ZfcLi 1 ( afc ) 2 ; we measure not 
the single operator n ■ T of the ancilla, but the N — 1 
operators specified by properly chosen unit vectors, cor- 
responding to N diagonal density matrices 



V, 



w = t 



2N 
N-l 



rpo 



(A.18) 



(with no summation over i on the right hand side), one 
obtains 



A' 



N 



2_/1 a Ip 2(N - 1) ^ 



\ ^ rpOirp{3 



N 



a u - N _ 1 



S^. (A.19) 
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If one prepares the ancilla in level i, measures its pop- 
(A 14) ulation after the gate operations and then sums up the 
results for all the measurements from i = 1 to N, the 
term quadratic in a reads (with summations made ex- 
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plicit here) 



level in p, Tr{(rii • T)p}, as 



V 



y^(g x m) ■ (a x m) 

L 

= x »»„) • (a x hp) 



i=l 



N 



8{N 



-^(a-r)^(a. 



t / /// 



w ( rp(y.rr\OL i rriQ; ryna ryrpOtrpQ \ 

x \ 1 U 1 U -r 1 mm 1 mm ~ A1 U 1 mm) 



N 



2{N 

N 
N - 



-{TT{(a-T) 2 }-J2(a-T) u (a-T) u 



a fc ) 2 = 



TV- 



(A.20) 



It is evident that together with the information on a Q , 
that can be obtained by measuring the population of each 



i=l 



1 * 

-^^Trj^.TH, 



(A.21) 



i=l 



one is able to estimate a 2 and thus purity II via (|A.8j) 
from the experimental data Tr{(n^ -T)a/ 2 )} and Tr{(ni • 
T)p} (i = 1, . . . , AT) (the populations of each level in u/ 2 ) 
and p). 

It is interesting to note that for an A^-level system pu- 
rity is expressed function of invariants (Casimirs) 
of SU(A r ). This is a peculiarity of the method proposed 
(VSWAP gates in succession, as in Fig. [T|), but might 
be of more general significance. We leave this issue for 
future investigation. 
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